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Abstract 

The electrical conduction properties of G4-DNA are investigated using a hybrid approach, which 
combines electronic structure calculations, molecular dynamics (MD) simulations and the formu- 
lation of an effective tight-binding model Hamiltonian. Charge transport is studied by computing 
transmission functions along the MD trajectories. Though G4-DNA is structurally more stable 
than dsDNA, our results strongly suggest that the potential improvement of the electrical trans- 
port properties in the former is not necessarily related to an increased stability, but rather to 
the fact that G4 is able to explore in its conformational space a larger number of charge-transfer 
active conformations. This in turn is a result of the non- negligible inter-strand matrix elements, 
which allow for additional charge transport pathways. The higher structural stability of G4 can 
however play an important role once the molecules are contacted by electrodes. In this case, G4 
may experience weaker structural distortions than dsDNA and thus preserve to a higher degree its 
conduction properties. 
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I. INTRODUCTION 



For the past two decades there has been a revival of interest in issues related to charge 
migration in DNA based oligomers. On one side, charge transfer (CT) in DNA is supposed 
to play a key role in self-repair processes of DNA damage in natural conditions, i.e. via 
oxidative stress pQ. On the other side, DNA may have a huge potential application field 
in the development of complex nano-scale electronic devices with self-assembling properties, 
since self-assembly, recognition and fully-automatic synthesis of oligonucleotides [2HD can 
allow for building almost any imaginable two and three dimensional DNA motifs (e.g. DNA- 
Origami) [5_H7]. However, besides considerable progress in elucidating the possible charge 
transfer mechanisms efficient in DNA, there is still disagreement about the conduction prop- 
erties of dsDNA, for several experiments obtained quite controversial results ranging from 
insulating [8], over semi-conducting [9] to even metallic-like behavior [101412] . This broad 
variability is related not only to the poor control of the experimental conditions in earlier 
transport experiments, but is also due to the intrinsic structural and electronic complex- 
ity of DNA oligomers (e.g. specific base pair sequence), which makes the formulation of a 
general theoretical framework very challenging. Nevertheless, several experiments appeared 
recently pT2HT¥j which have shown electrical currents in the nA range despite the differences 
in the studied base sequences and in the experimental setups. One central issue which has 
emerged by the theoretical treatment of charge transport is the necessity to consider the 
conformational fluctuations of the molecular frame not as a weak perturbation but rather 
as a crucial factor promoting or hindering charge motion [15-20] . 

It has been meanwhile shown that it is possible to synthesize DNA derivatives, which do 
not necessarily have the double-strand structure of natural DNA. Thus, C and G-rich DNA 
strands are able to form four-stranded quadruplex structures, like the i-motif structure 
composed of two parallel hemiprotonated duplexes intercalated into each other into a head- 
to-tail orientation, and the G4-quadruplex (also known as G4-DNA) which is formed by 
either one, two or four G-rich DNA strands in a parallel or anti-parallel orientation [2T] . 
The latter is supposed to play an important role in some biological processes such as in 
telomeric DNA regions, where it inhibits telomerase and HIV integrase [22]. Additionally, 
G4-DNA is known to interact with various cell proteins that cause diseases like Bloom's and 
Werner's syndromes [23] • Moreover, G4-DNA is cytotoxic towards tumor cells and therefore, 
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might be a key for the design of anti-cancer drugs [21]. Eventually, G-quadruplexes are found 
to be thermally more stable than dsDNA [25J. Gilbert and Sen resolved the x-ray structure 
of the first G-quadruplex and proposed the formation of these unique structures in presence 
of monovalent alkali ions [26H28] . see also Fig. [l](a). Consequently, G4-DNA can be regarded 
as stacks composed of individual planar G-tetrads (or G-quartets) as shown in Fig.fiTb), each 
formed by 4 guanines bound together by 8 hydrogen bonds via Hoogsteen pairing [29]. For 
detailed reviews on the structural diversity and special mechanical properties of G4-DNA 
based molecules see Ref . [301 EI] • 

Recently, Porath et al. [22] have shown the larger polarizability of G4-DNA compared to 
dsDNA, when attached to gold surfaces. This might be an indicator of improved CT prop- 
erties and the potential advantage of using G4-DNA rather than dsDNA in molecular scale 
devices. Eventually, the higher conductance in G4-DNA could be attributed to the increased 
structural stability and higher number of overlapping 7r-orbitals. Furthermore, experimental 
work by Kotlyar et al. [331 El] revealed that up to 300 nm long G4- wires could be synthesized 
revealing a considerable stability even in the absence of metal ions. 

Also, theoretical studies have been performed on the stability and rigidity of G4-DNA. 
Molecular dynamics (MD) simulations by Spackova et al. [351 EE] could confirm the experi- 
mental results. For one thing, G4-DNA is more rigid than dsDNA and for another monova- 
lent ions within the quadruplex cavity are necessary for the stability of short G4-DNA. On 
the other hand, MD simulations [37] pointed out the stability of longer G4 molecules (with 
24 G-tetrads) in absence of metal ions which is consistent with the synthetic procedure of 
ref [331 El]. Electronic structure calculations by Di Felice and co-workers [38H13] showed a 
higher degree of derealization of the electronic states in G4 than in dsDNA. Moreover, the 
presence of metal ions may contribute additional states supporting CT, though this is a still 
unresolved issue. 

Guo et al. [lH US] have recently shown using the Landauer theory (coherent transport) that 
G4-DNA exhibits much larger derealization lengths at the band center compared to double- 
stranded poly(G) DNA. Surprisingly, it was also found out that the derealization length 
can be even enhanced via environment-induced disorder through the backbones. Though 
these studies assume a static atomic structure, they nevertheless suggest that disorder plays 
an important, non-trivial role in mediating CT in DNA. 

Similarly, it is known from other biomolecules such as proteins that electron transport is 
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dominated by non-equilibrium fluctuations which has lately been analyzed by Balabin et 
al. [IB] . As demonstrated in previous studies, idealized static structures are not representa- 
tive when considering CT properties, since dynamical as well as environmental effects were 
shown to be too important [HH5T]- It was indicated that dynamic disorder has dramatic 
effects, since it suppresses CT in homogeneous sequences on the one hand, but can enhance 
CT in heterogeneous sequences, on the other. For instance, in our previous work [20J, the 
conductance of the Dickerson dodecamer (sequence: 5'-CGCGAATTCGCG-3') was found 
to be almost one order of magnitude larger in solution (QM/MM) than in vacuo. Interest- 
ingly, similar findings have recently been shown by Scheer and co-workers, for they obtained 
for a heterogeneous sequence (i.e. A and G bases are present) with 31 base pairs [68] a 
two order of magnitude larger conductance in solution compared to the in vacuo measure- 
ments [13]. Furthermore, it could be shown that only a minor part of the conformations is 
CT-active [20]. Therefore, neglecting these significant factors or assuming a purely random 
disorder distribution can lead to a considerable loss of a vital part of CT-relevant structural 
and electronic information. 

Relying on our previously developed methodologies to deal with CT in different DNA 
oligomers [HI [20l [52] , we present here a detailed investigation of charge transport in G4- 
DNA. 

The paper is organized as follows, in Sec. [TT] details of the MD simulations and of the elec- 
tronic structure approach are described. Subsequently, the MD and the electronic structure 



data are analyzed in Sec. Ill A and III B , respectively. Finally in Sec. Ill D , CT results for G4 
molecules are compared with those of dsDNA in solution and in vacuo. More importantly, 
we point out significant factors which are responsible for the enhanced conductance in G4 
molecular wires. 



II. METHODOLOGY 

A. Starting Structures and Simulation Set-Up 

The molecules used in this work are based on the x-ray crystal structure of a tetrameric 
parallel-stranded quadruplex (244d) formed by the hexanucleotide sequence d(TG4T) in the 
presence of sodium ions. The structure has been resolved at 1.2 A resolution by Laughlan 
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et al. [53]. As shown in Fig. [j](a), it contains 2 pairs of quadruplexes (TG4T) 4 , for which 
each pair is stacked coaxially with opposite polarity at the 5' ends. Moreover, 9 sodium ions 
are located inside each stack of quadruplexes illustrating well ordered G4-DNA constructs. 
However, the terminal thymine residues, shown in gray, were not completely resolved because 
of high thermal disorder [53]. As indicated in Fig. [I] only one of these 4 quadruplexes is 
taken as our basis structure, which corresponds to the parallel strands A, B, C and D 
from the PDB file 244d. This quadruplex, from now on denoted as (TG4T) 4 , will be used 
for simulations and calculations. Furthermore, two longer G4 quadruplexes with 12 and 
30 tetrads, denoted as (Gi 2 ) 4 and (G 30 )4, are generated by omitting the terminal thymine 
residues and adding subsequently G4 tetrads with a distance of 3.4 A and twisted by 30°. 
For comparison, corresponding double-stranded B-DNAs with base sequence 5'-TGGGGT- 
3', poly(G) denoted as G12 and G30, respectively, and a sequence containing 31 base pairs [68] 
(Scheer) are built with the make-na server [54] . 

Depending on the helix length the G4 and dsDNA molecules are solvated in a rectangular 
box with 4000 to 11000 water molecules using the TIP3P model [55] and periodic boundary 
conditions (PBC) are applied. In order to neutralize the total charge of the system, due to 
the negatively charged backbones, an appropriate number of sodium counterions (Na+) are 
added. Concerning the quadruplexes, we carried out 4 simulations for (TG4T) 4 and (612)4, 
respectively, firstly without any central monovalent alkali ions within the quadruplex, and 
then in presence of either lithium (Li+), sodium (Na+), or potassium (K+) ions within 
the respective G4 molecules. Therefore, either 3 ((TG 4 T) 4 ) or 11 ((Gi 2 ) 4 ) of these ions 
are subsequently placed in the center between two G4 tetrads coordinated by 06 oxygen 
atoms (see also Fig. [jjb)). Nevertheless, the longer (630)4 quadruplex is simulated only 
twice in absence and presence of central sodium ions. All simulations are carried out with 
the GROMACS software package [56J using the AMBER parm99 forcefield [57] including 
the parmBSCO extension [SB]. After a standard heating-minimization protocol followed by 
a 1 ns equilibration phase, which is discarded afterwards, we performed 30 ns (50 ns for 
(630)4) MD simulations with a time step of 2 fs. Snapshots of the molecular structures 
were saved every 1 ps, for which the charge transfer parameters were calculated with the 



SCC-DFTB-FO approach as described in Sec. II B 
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FIG. 1: a) Tetrameric G4 x-ray crystal structure (244d [53]) with four units (TG4T)4. Backbones 
are indicated as light blue ribbons, central coordinated sodium ions as dark blue spheres, terminal 
thymine and guanine residues are shown in gray and orange, respectively. The highlighted quadru- 
plex is used as starting structure and also to generate longer G4 molecules (612)4 and (630)4. b) 
Lewis structure of a single G4-tetrad with C^h symmetry containing a monovalent metal ion in its 
center. Guanines are bound together by 8 hydrogen bonds via Hoogsteen pairing |29j . The metal 
ion is coordinated either coplanar by 4 or cubic by 8 06 oxygen atoms of the respective guanines 
which depends on the ionic radius of M + . 

B. Electronic structure 

In this section we will shortly describe how the electronic structure of G4-DNA molecules is 
mapped to a coarse-grained transfer Hamiltonian using the fragment orbital (FO) approach. 
The method has already successfully been applied to dsDNA molecules [191 EH E2J [59] . A 
more detailed description of the methodology can be found in Ref. [60]. To begin with, 
the electronic structure of the DNA system is written in an effective tight-binding basis 
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(fragment basis) as: 



H = '^^e i ala i + '^^T i j(a\aj + h.c.) 



(1) 



The onsite energies £j and the nearest-neighbor hopping integrals characterize, respec- 
tively, effective ionization energies and electronic couplings of the molecular fragments . The 
evaluation of these parameters can be done very efficiently using the SCC-DFTB method [61] 
combined with a fragment orbital (FO) approach [17]: 



and 



Si = - 



ij 



H, 



KS 



H: 



KS 



(2) 



(3) 



The indexes i and j correspond to the fragments of the quadruplexes which, in this model, 
are constituted by single guanine bases. Accordingly, the molecular orbitals fa an d fa are 
the respective highest-occupied molecular orbitals which are obtained by performing SCC- 
DFTB calculations for these isolated fragments, i.e. the individual guanine bases. Using a 
linear combinations of atomic orbitals ansatz <pi = Y2 c^rj^ the coupling and overlap integrals 
can be efficiently evaluated as: 



H 



KS 



fJ,U 



and 



Sij — c^cl (r]fj,\rj v ) — c^c^S^u. 



(4) 



(5) 
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H^ v and are the Hamilton and overlap matrices in the atomic basis set as evaluated 
with the SCC-DFTB method. Note, since is built from non-orthogonal orbitals fa and 
(f)j, which is for various problems not suitable, we apply the Loewdin transformation [62] : 



S-fT S _ t 



(6) 



The effect of the environment, i.e. the electrostatic field of the DNA backbone, the water 
molecules and the counter-ions (including the central alkali ions within the quadruplex), is 
taken into account through the following QM/MM Hamiltonian: 

1 



H, 



flU 



H 



Aq s (<y a5 + 7/35) 




(7) 



Aqs are the Mulliken charges in the QM region and Qa are the charges in the MM region, 
i.e. the DNA backbone, counterions and water molecules. The coupling to the environment 
is therefore explicitly described via the interactions with the Qa charges. In the following, 
we will denote the calculation set-up based on the complete expression in Eq. [T]as QM/MM; 
neglecting the last term will be denoted as "vacuo" . 

The electronic parameters e, and T% are evaluated for every snapshot of the simulations. 
Subsequently, the transfer Hamiltonian is used to calculate the CT in the quadruplexes as 
described in Sec. Ill CI 

C. Charge transport through a 4-stranded quadruplex 

Once the transfer Hamiltonian has been constructed, the transport properties will be cal- 
culated using Landauer theory for each snapshot along the MD trajectories. For this, we 
consider a two-terminal set up where the G4 oligomer is contacted to left (L) and right (R) 
metallic electrodes through the four terminal bases. A central quantity to be computed is 
the molecular (G4) Green function, which can be obtained through a Dyson equation: 

G- 1 ^) = £1-H-£ L -£ R . (8) 

Here El and are self-energy matrices characterizing the coupling to the electrodes. To 
simplify the calculations we do not take the full energy dependence of the self-energies into 
account, but rather use the so called wide band limit, where T, L and T, R are replaced by 
energy-independent parameters [Ml El] : 

(S L )y = -iTL^ifc^-fc, {k = 1,2,3,4) 

(E R ) y = -ij R 6 lh 6 jk ,(k = N-3,N-2,N-l,N). 

We set 7l and 7r to 1 meV. Within the Landauer approach, the transmission function T(E) 
for a given set of electronic parameters is then obtained as: 

T(E) = Tr [r L Gr L G+] (9) 

where Tl and Tl are the broadening matrices calculated as the anti-hermitian part of El/r : 

(10) 



L/R 



S L /R - S+ /R 
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Using the former expressions, the conformational (time) dependent electrical current will be 
simply calculated by: 

- /(s-S,+ y) T(S,t) (11) 

The I-U characteristics presented in this work should be interpreted only qualitatively, for 
the Fermi energy (Ep) is artificially placed as average of the onsite energies for all guanine 
sites and for each snapshot, respectively. The reader should note that the current could 
exhibit quite different shapes depending on where Ep is located. However, only the current- 
voltage gap is affected by Ep, whereas the maximum current obtained at high voltages is 
not altered. 

The transmission function and the current are evaluated for every snapshot of the MD 
simulation in order to obtain statistical average quantities. However, we are aware of the 
limits of the coherent transport model which is valid only in the adiabatic regime. Hence, we 
assume that the timescales of the CT process is shorter than the fastest dominant structural 
fluctuations. This issue has also been addressed in Ref. [20J. 



III. RESULTS 



A. Structural stability and rigidity in parallel stranded quadruplexes 

First of all, the (TG4T)4 quadruplex from the crystal structure 244d has been simulated in 
order to validate further simulations with 12 and 30 tetrads. Corresponding RMSD time 
series in absence and presence of central ions as well as molecular snapshots after 30 ns 
simulation time can be found in Sec. SLA. in Ref. [72]. Based on the crystal structure 244d 
we built up an artificial all-parallel stranded guanine quadruplex (612)4 without terminal 



thymine residues and 12 stacked guanine tetrads (for details see II A). As for (TG4T)4, we 
carried out 4 simulations, one in absence and three in presence of 11 lithium, sodium and 
potassium ions initially placed in the center of two adjacent G4 tetrads. The RMSD time 
series and molecular snapshots after 30 ns MD simulation for (612)4 and poly(G) (G12) are 
given in Fig. [2] and Fig. |3j respectively. 
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FIG. 2: RMSD of G4 and poly(G) DNA: a) parallel stranded quadruplex (Gi2)4 in absence and 
presence of centrals ions Li + , Na + and K + , b) comparison with double-stranded Gi2- 

The RMSD values show that without central ions the quadruplex cannot be considered 
as a G4 structure in equilibrium. In fact, the RMSD increases for the whole simulation 
time from 1.5 to 3 A indicating considerable structural changes. In contrast, the RMSD 
for structures containing central ions reveal rather consistent values around 1 to 1.4 A with 
only minor fluctuations suggesting rigid G4 structures with only slight structural deviations. 
Compared to the (TG4T) 4 simulations, the type of ions within the quadruplex has a lesser 
influence on the structural deviations, since the RMSD values are rather similar. In order to 
validate these small structural deviations and therefore the rigidity of G4-DNA we carried 
out a RMSD calculation for double-stranded Gi 2 (Fig. [21(b)) . Generally, the RMSD values 
and their fluctuations for the dsDNA molecule are substantially larger than those for the 
quadruplex structures. This clearly supports the notion that guanine quadruplexes are 
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FIG. 3: Molecular snapshot of parallel stranded quadruplex (Gi2)4 after 30 ns MD simulation in 
absence and presence of centrals ions Li + (green spheres), Na + (blue spheres) and K + (purple 
spheres). 

much more rigid than corresponding dsDNA and therefore G4-DNA reveals less dynamical 
disorder which might promote CT. 

The snapshots in Fig. [3] exhibit that (612)4 m absence of central ions is not completely 
destabilized which has not been the case for (TG4T)4 (see Sec. SLA. in Ref. [72J). The inner 
tetrads roughly maintain a regular G4 structure, whereas the outer ones show a considerable 
disordered behavior with no planarity at all, even the G4 hydrogen bridge pattern has been 
entirely destroyed. Recently, Cavallari et al. have shown in their MD simulations that long 
G4 quadruplexes can be stable without central ions |37]. This finding was supported by 
experimental results in which longer G4 quadruplexes could be synthesized and characterized 
in absence of monovalent ions 1 33 1 . 




Therefore, an additional 50 ns simulation of a quadruplex with 30 tetrads (630)4 in absence 
of central ions is carried out. The simulations clearly exhibit that the major part of the 
quadruplex forms a quite regular structure with adequate parallel stacking. Therefore, it 
seems to confirm previous results that large G4 wires in absence of central ions become more 
stable, reasonably due to increased stacking interaction. However, like for the dodecamer, 
towards the ends the quadruplex is slightly deformed, especially at the 3' terminus. More- 
over, RMSD values reveal that the quadruplex is not completely structural balanced. The 
major part of structural fluctuation originates from the 5' and 3' ends which is shown by four 
groups of large peaks in root mean square fluctuations per atom (RMSF). The corresponding 
data is shown in Fig. S2 in Ref. [72]. 

The quadruplexes which contain central monovalent ions in Fig. [3] exhibit highly regular G4 
structures indicating again the strong rigidity and stability of these molecules. Generally, 
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the ions appear to prefer different coordinations within the quadruplex cavity, that is lithium 
favors the planar position within the tetrads while sodium and potassium are most likely 
to be found in the center of two G4 planes. However, one sodium ion in the middle of the 
quadruplex is located in-plane and two others are rather closer to one of the corresponding 
tetrads than sandwiched between them. Also, lithium ions are sometimes found slightly 
apart from the in-plane position. This might be an indicator for the longitudinal mobility of 
central sodium and lithium ions within guanine quadruplexes. In contrast, potassium ions 
are rather fixed in their center position between two tetrads, most likely because of their 
larger ionic radii compared with lithium and sodium. The resulting larger rigidity of G4 
compared to dsDNA plus the preferred locations of central ions within the quadruplex are 
in perfect agreement with findings obtained by Spackova et al. as well as previous work by 
Cavallari et al. [351 EH [65]. 

Notwithstanding, the ions are equally capable of moving out of the quadruplex through the 
minor groove and the 5' and 3' ends. In our simulations 3 sodium, 3 lithium and 2 potassium 
ions went out of the G4 dodecamer, respectively. Sometimes, these ions remain in the minor 
groove for several ps. Likewise, solvent molecules are able to infiltrate the quadruplex via the 
minor groove and the 5' and 3' ends provided that these locations are not already occupied 
by initially placed ions. Usually, the infiltration of solvent molecules into the helix will 
not significantly destabilize the structure as long as there is at least a minimum amount of 
central monovalent cations left in the quadruplex. Otherwise, the interior solvation might 
destroy the regular 4-stranded structure as can be seen for the structures without central 
ions in Fig. [3j These observations are consistent with results of previous works [351 [33, [65] . 

B. Electronic Structure of G4 Tetrads 

1. Idealized Model 

Before calculating CT characteristics, we have to focus on the electronic parameters, espe- 
cially the molecular orbitals involved in the hole transfer process. Similarly to our previous 
work on CT parameters in dsDNA [HO] we carried out benchmark calculation of G4 tetrads 
with DFTB and compared to DFT (PBE and B3LYP) as well as HF calculation. Since we 
intend to deal with hole transfer we are most notably interested in the highest molecular 
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TABLE I: Energies of highest occupied molecular orbitals for a idealized G4 tetrad with C^- 
symmetry, comparison between DFTB, DFT and HF, for the latter the 6-31G(d,p) basis set is 
used, all values in eV. 



Method HOMO HOMO-1 HOMO-2 HOMO-3 HOMO-4 



DFTB 


-4.588 


-4.593 


-4.593 


-4.597 


-4.870 


PBE 


-4.356 


-4.383 


-4.383 


-4.408 


-5.267 


B3LYP 


-5.170 


-5.196 


-5.196 


-5.222 


-6.637 


HF 


-7.777 


-7.803 


-7.803 


-7.828 


-10.827 



orbitals which ought to have 7r-symmetry in order to have sufficient MO overlap along the 
quadruplex. Tab. [T] shows energies of the 5 highest occupied MOs for a idealized G4 struc- 
ture with C^-symmetry obtained after geometry optimization with B3LYP/6-31G(d,p). 
Corresponding visualization of MOs can be found in Tab. S3 in Ref. [72] . 
The first four occupied MOs, HOMO to HOMO-3 are very close in energy for all methods 
respectively and they have 7r-symmetry. The range is about 0.01 eV for DFTB and 0.05 
eV for the rest. In contrast the HOMO-4 is separated significantly by a gap of about 
0.3 eV for DFTB and at least more than 0.8 eV for the other methods. For DFT based 
methods this MO has cr-symmetry. As a result HOMO to HOMO-3 are the orbitals used for 
the transport calculations, for the close energy range between them supports the notion of 
band-like transport provided large G4 stacks are considered. By contrast, static disorder in 
dsDNA due to differences in ionization potentials between four different nucleotides (A,C,G 
and T) leads to large energy gaps for tunneling and therefore to rather localized states. 
Basically DFTB MO energies lie between those for PBE and B3LYP. 

More importantly, these MOs appear to be linear combinations of HOMOs for the isolated 
guanine molecules which can be shown by comparison with our previous results [60J . There- 
fore, we are able to map the electronic structure onto single guanines rather than onto whole 
tetrads, since this reduces the computational costs immensely. 

In Sec. II. B. in Ref. [72], we discuss the effect of central ions on the MOs of G4-tetrads both 
purely electrostatically and quantum mechanically. As a main result, we find that the MO's 
are shifted down due to the electrostatic interactions with the ions. This effect is captured 
by a classical treatment of the ions, therefore, in our MD simulations we can treat the ions 
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TABLE II: Energies of highest occupied molecular orbitals of a single tetrad of (TG4T)4 for a 
MD snapshot containing the full electrostatic environment built of negatively charged backbone, 
solvent and counterions. Comparison between DFTB, DFT and HF, for the latter the 6-31G(d,p) 
basis set is used, all values in eV. 



Method HOMO HOMO-1 HOMO-2 HOMO-3 HOMO-4 



DFTB 


-5.390 


-5.409 


-5.467 


-5.544 


-5.756 


PBE 


-5.052 


-5.156 


-5.206 


-5.273 


-5.986 


B3LYP 


-5.871 


-5.949 


-5.994 


-6.121 


-7.356 


HF 


-8.322 


-8.328 


-8.362 


-8.711 


-11.113 



with MM rather than with QM. 
2. Dynamics and QM/MM 

Up to now we only calculated ideal G4 tetrads which where highly symmetric, therefore we 
applied the previous routine on a molecular snapshot taken from a MD simulation. That 
way we can analyze the effects of dynamics as well as the purely electrostatic environment 
built of backbone and solvent charges on the electronic structure of G4 tetrads. As shown in 
Sec. II. B in Ref. [72] single alkali ions placed in the quadruplex cavity possess a substantial 
impact on the MO energies of hole transfer states. Therefore, it is vital to include the MM 
environment in our calculations. The MO energies in Tab. [IT] for all applied methods are 
shifted down by about 0.5 to 0.8 eV compared to those for the idealized tetrad without MM 
environment. Also, the energy range in which HOMO to HOMO-3 are located has become 
considerably larger from 0.01 to about 0.15 eV for DFTB and equally from 0.05 to almost 
0.4 eV for DFT and HF. Accordingly, these hole transfer states become rather localized onto 
the single guanines than onto the whole tetrads as soon as the electrostatic environment is 
considered, although the energy range is still very small. 

In Tab. S5 in Ref. [72] we provide snapshots of corresponding MOs for these states. However, 
this effect of localization clearly supports the notion to build up a transport model with single 
guanines as fragments rather than whole tetrads. It has not only lower computational costs, 
for only 16 instead of 64 atoms have to be calculated per fragment, but it also offers a large 
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FIG. 4: Scheme for electronic couplings in G4-DNA shown for the two inner most tetrads in 
(TG4T)4. Representatively, Tl indicates intrastrand, T2 interstrand and T3 in-plane couplings. 

flexibility because CT can occur through a multitude of pathways along the quadruplex. 



C. Electronic Parameters 



In this section we analyze CT parameters obtained for the simulations of (TG4T)4 and 
compare them to those for the crystal structure 244D as well as for an idealized G4 stack. 



The parameters were calculated as described in Sec. II B Before showing the results, first the 
applied fragment methodology is introduced. The scheme in Fig. [4] shows different types of 
couplings present in G4-DNA. Here T3 represents electronic couplings within the respective 
G4-tetrads in-plane. Whereas Tl and T2 denote intra- and interstrand couplings lengthways 
to the quadruplexes which occur on either one strand or between two strands, respectively. 
In contrast to dsDNA, experiments suggest that also competing horizontal CT can occur in 
G4-DNA jSg. 

Clearly, a charge can follow several pathways along the quadruplex strongly dependent on 
those three couplings. This should be an advantage compared to dsDNA. Even in small 
G4 stacks there is a large number of electronic couplings. We will see that only a minor 
part of them will be vital for CT in G4-DNA. For instance, T3 couplings are usually quite 
small, especially the diagonal in-plane couplings in Fig. [5] (G1-G3,G2-G4, ...) are generally 
negligible. Also, most of T2 couplings are small, except those between adjacent strands (Gl- 
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TABLE III: Average onsite energies (e) and electronic couplings Tl, T2 and T3 with standard 
deviations for the two central guanine tetrads of (TG4T)4. Energies obtained from MD simulations 
in absence and presence of ions are compared to the crystal structure (244D), an idealized G4 
dimer stack, and also to the corresponding dsDNA structure, i.e. the two central guanines in 5'- 
TGGGGT-3'. Note that averaging is carried out not only along the MD time series but also over 
8 and 2 G bases for G4 and dsDNA, respectively. All values in eV. 
Type ideal 244D MD no ions MD Li+ MD Na+ MD K+ dsDNA 
e -4.895 -4.905±0.062 -4.812±0.368 -5.339±0.370 -5.400±0.350 -5.201±0.354 -4.790±0.371 
Tl 0.028 0.051±0.011 0.039±0.028 0.031±0.021 0.031±0.021 0.029±0.020 0.052±0.034 
T2 0.001 0.012±0.002 0.010±0.013 0.022±0.014 0.015±0.012 0.013±0.010 0.004±0.005 a 
T3 0.009 0.007±0.003 0.009±0.009 0.007±0.005 0.006±0.004 0.007±0.004 0.012±0.008 b 

a interstrand coupling for G\C, G/C=0.013±0.014 eV 

h Coupling within the WCP between G and C, note there is an energy gap of 0.4 eV 

G6,G2-G7, ...), since these are rather close to each other. The largest couplings are certainly 
found for Tl (intrastrand) and they are well comparable to those for double-stranded poly(G) 
DNA [HO]. 

A summary of the average electronic couplings and onsite energies for the simulations of 
(TG4T)4 compared to the crystal structure (244D) as well as to the corresponding dsDNA, 



i.e. the two central guanines in 5'-TGGGGT-3' is given in Tab. [ITT 

As has been discussed above, the central ions stabilize the G4 structure significantly, the 
structural fluctuations of G4 with ions are much lower compared to G4 without ions or 
dsDNA. The effect of ions on the electronic parameters is twofold: First of all, the onsite 
energies are shifted down by about 0.5 eV but , surprisingly, the enhanced structural stability 
does not lead to smaller fluctuations of the onsite parameters. They are in the order of 0.4 
eV, similar to the situations in dsDNA, as discussed recently (59]. These large parameter 
variations are introduced by solvent fluctuations (and not by fluctuations of the DNA/G4 
structure itself), therefore, they are not affected by the higher structural rigidity. The Tl 
values, on the other hand, are even lower in the G4 structures with ions. This would indicate, 
that these structures conduct even less, when compared to four strands of dsDNA. However, 
the T2 and T3 values in G4 are still large, indicating that interstrand transfer can occur 
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quite frequently. This opens a multitude of pathways for charge transport in G4 (compared 
to dsDNA), which will be the key to understand G4 conductivity, as discussed in more detail 
below. 

As a first result, we do not see any indication that the higher structural stability of G4 
with central ions will lead to a higher conductivity due to a reduced dynamical onsite 
disorder or due to increased electronic couplings, i.e. because of somehow better stacking 
interactions owing to the more regular structure. Therefore, the higher conductivity must 
have different reasons. The more stable structure leads to smaller couplings, in contrast to 
prior expectations [22] • A more detailed analysis of onsite energies and electronic couplings 
is given in Tab. S6-S11 in Ref. [72J. 



D. Coherent Transport in G4-DNA 

The parameters e and Ty analyzed in the previous section are now used to evaluate the 



charge transport properties using the approach described in Sec. |II C[ First, we focus on 
static structures to get insight in particular CT features for G4 quadruplexes and also to have 
a reference for the following calculations which are applied to the whole MD trajectory. Later 
on, we are interested in ensemble averages, for single snapshots or static structure cannot 
elucidate the CT process in DNA. Subsequently, we will compare the CT properties of G4- 
DNA with those of double-stranded DNA and distinguish the differences in CT efficiency 
by means of conformational analysis. Additionally, we determine the effects of the MM 
environment and in particular of the central ions on CT in G4 quadruplexes. 



1. Static Structures 



Primarily, we compare T(E) and I(U) for three static structures: i) idealized G4-Dimer as 



used in Tab. Ill, ii) the two central tetrads of the crystal structure 244D and iii) double- 
stranded idealized B-DNA (G2), i.e. two stacked G-C base pairs (36° twist and 3.4 A 
rise). As can be seen from Fig. [5] there are only slight structural differences between the 
idealized G4-dimer and the stacked tetrads of the crystal structure 244D. Fig. [6] exhibits 
the corresponding transmission functions showing eight resonances due to the eight G sites 
of the transport model. However, since the idealized G4-dimer has a very regular structure 
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FIG. 5: Static structures, left: Idealized G4-dimer (ideal), 30° twist and 3.4 A rise; middle: G4 
tetrads of crystal structure 244D (244d); right: double-stranded idealized B-DNA (G2), i.e. two 
stacked G-C base pairs (36° twist and 3.4 A rise). 

these eight peaks lie in a very narrow energy range of about 0.15 eV and have similar high 
transmission values of at least 0.9. On the other hand, the eight peaks of the two central 
tetrads of 244D lie in a broader energy range of about 0.26 eV and have lower transmission 
values ranging from 0.4 to 0.8. This indicates, how small structural deviations lead to a 
broadening and a decrease of the transmission function. Nevertheless, we know from the 



values in Tab. |III| that the Tl and T2 couplings are much larger for the crystal structure 
than for the idealized structure, thus the cause for decreased conductivity has to be the 
onsite energy disorder. 

Fig. [6] also shows the corresponding current- volt age characteristics. The current for the G4 
structures is significantly larger than for G2. This is what we would expect, because we 
compare two idealized molecular units, one with two relevant transmissive states and the 
other with eight within the relevant energy window. Note that in all current calculations 
shown in this work, the Fermi energy is artificially placed as the mean value of onsite energies 
for each snapshot respectively. 

Generally, the conductance for the idealized double-stranded poly(G) DNA would become 
twice as large once the cytosine states are included, though it is well established that the 
C states have only a minor impact in hole transfer due to their higher oxidation potential 
than G. 



2. Influence of Dynamics and Environment 

Previous studies revealed that idealized static structures cannot exhibit reasonable CT prop- 
erties in DNA, since dynamical as well as environmental effects were shown to be crucial 
[4~TH5"T] . Dramatic effects are induced by dynamic disorder, for it suppresses CT in homo- 
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FIG. 6: Charge transport in static structures from Fig. [5j comparison of transmission function 
T(E) and current I(U) for idealized G4-dimer, the two central tetrads of crystal structure 244D 
and poly(G) (G2). Note, only the stacked G bases are considered, whereas the cytosine bases and 
the backbone composed of sugar and phosphate groups are neglected. 

geneous sequences on the one hand, but can also enhance CT in heterogeneous (random) 
sequences on the other. Moreover, only a minority of conformations appears to be CT-active 
as has been indicated in Ref. [20]. Therefore, neglecting these significant factors or assuming 
purely random distributions for dynamical disorder leads to a considerable loss of a vital 
part of the CT in DNA. On this account the CT properties are evaluated for every snapshot 
along classical MD trajectories which then leads to ensemble averaged quantities, that is the 
average transmission function (T(E)) and the current (I(U)). 

In Fig. [7] (T(E)) and (I(U)) are shown for both quadruplex molecules (with central sodium 
ions), the two central tetrads of (TG4T) 4 (a) and (G§)4 (b) as well as their corresponding 
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double-stranded poly(G) analogues G2 and Gg[69j. Because of the substantial onsite energy 
fluctuations of about 0.4 eV due to the DNA environment and dynamics, the transmission 
function reveals large broadening for both DNA species. The transmission maxima for 
the quadruplexes are shifted to lower energies by about 0.3 eV due to the presence of the 
central sodium ions. Basically, the average transmission is strongly reduced compared to 
the idealized static structures in Fig. [6j Nevertheless, the maximum of (T(E)) for the two 
central tetrads in (TG4T) 4 is almost five times larger compared to G2. Accordingly, the 
maximum current is about 4.4 times larger in the quadruplex. 

Considering the octamers in Fig. ^b) the conductance difference between G4 and poly(G) 
even increases. Here, the average transmission in the relevant energy range for the quadru- 
plex is to a great extent two orders of magnitude larger than those for the poly(G) sequence. 
The poly(G) spectra shows much larger spikes at certain energies. These spikes indicate the 
strong impact of dynamics and may be explained by the existence of few charge transfer ac- 
tive conformations, which dominate the average transmission, i.e. few conformations, which 
feature a high transmission. This may indicate that even longer sampling (more than 30 ns) 
would be required to converge the spectra. Notwithstanding, it also reveals that in dsDNA 
the average current is dominated to a much larger degree on few non-equilibrium structures, 
as concluded from our earlier work [20]. As a result, the average current (I(U)) for (Gs)a is 
almost two orders of magnitude larger than for Gg suggesting that the enhancement of CT 
in G4 with respect to poly(G) might grow with increasing DNA length. 
Thus, additional sets of CT calculations for longer DNA species with 14 and 20 tetrads 
(base pairs), respectively, are carried out. The corresponding data was obtained from MD 
simulations of a quadruplex (630)4 (including central Na + ), and double-stranded DNAs 
G30 and a heterogeneous sequence containing 31 sites [68] recently used by Scheer and co- 
workers in a CT-measurement [13] (Note, for CT calculations only the 14 (20) central sites 
are used). As expected, the transmission strongly decreases for both DNA species. However, 
in G4 this effect is not as strong as in dsDNA. For instance the current for (614)4 is only 
about 1 order of magnitude lower than for (Gs)4, although expanding to (620)4 the current 
drops significantly additional 3 orders of magnitude. In contrast, (I(U)) for poly(G) and the 
"Scheer" sequence decreases by more than 10 orders of magnitude by increasing the number 
of base pairs from 14 to 20. This indicate there is a much stronger distance dependence 
of CT in dsDNA, hence the notion of coherent CT for longer molecular wires might be 
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FIG. 7: Average transmission (T(E)) and current (I) obtained from MD simulation: Comparison 
between G-quadruplex and double-stranded poly(G) DNA. a) (G^i (Na + ) and G2, i.e. the two 
central tetrads (base pairs) of (TG4T)4 (5'-TGGGGT-3'), respectively, b) (Gg)4 (Na + ) and Gs- 
Generally, the two last tetrads (base pairs) at the 5' and '3 end are not considered for CT calcu- 
lations to avoid end effects, although the simulations were done with 12 tetrads and base pairs for 
G4 and dsDNA, respectively. 

considerably more likely in G4 than in dsDNA. Admittedly, the Landauer formalism used 
in this work performs well for short DNA species (less than 10 sites), where the transport 
is assumed to be at least partially coherent. On the other hand, it clearly fails for long 
DNA sequences, hence the CT results for the longer molecules should be interpreted only 
qualitatively and with caution. For instance, the currents obtained for the 14mer and 20mer 
of both dsDNA molecules poly(G) and the Scheer sequence are orders of magnitude smaller 



22 



than pA which is far beyond any measurable range. The complete data is given in Fig. S3 
in Ref. [72]. 



3. Analysis of CT Differences in G4 and dsDNA 

The significant conductance difference of G4 and dsDNA may not be attributed to the fact 
that G4 is composed of four poly(G) like wires. To analyze this further the CT in (G§)4 
with central sodium ions (full MD) is compared to various models, in which i) only intra- 
and interstrand[70j couplings are non-zero (intra+inter) and ii) the quadruplex is separated 
into its 4 single strands, for which the intrastrand transport is calculated independently 
and added up afterwards (S4s). Furthermore the results are analyzed with reference to G§ 
(4xpG). Note, for comparison the CT quantities of poly(G) are multiplied by 4. 
As appears from Fig. |8j the average transmissions for the two models: full MD and in- 
tra+inter in the spectral support region clearly reveal the largest plateaus, also showing 
similar peak structures. However, both quadruplex models exhibit similar moderate fluc- 
tuations. By contrast, the sum of the four single G4 strands (S4s) shows much larger 
fluctuations, comparable with those for G§ (4xpG). Moreover, the average transmission for 
E4s is significantly reduced and is even slightly lower than 4xpG. Interestingly, in S4s are 
barely CT-active conformations, i.e. single dominating peaks like in 4xpG which some- 
times even outreach the maximum transmission of the quadruplex (full MD). This might 
reflect the structural differences between quadruplexes and dsDNA, since the four strands 
in the quadruplex are not as flexible as those in poly(G) (see also RMSD fluctuations in 
Fig. [2]), hence in poly(G) the structural phase space is much larger and therefore, several 
high-transmissive structures can arise. On the other hand, the more rigid quadruplex covers 
only a smaller conformational phase space which indeed ensures a large number of structures 
showing moderate CT properties for each single strand, respectively. This clearly indicates, 
that the most important factor for the enhanced conductance in G4 are the interstrand cou- 
plings between the 4 strands in the quadruplex. Thus, if there are conformations in which 
the 4 isolated channels are not transmissive, there is a considerable probability that CT 



might occur via coupling between the individual stands. As can be extracted from Tab. Ill 
those interstrand couplings are sufficiently large with about 0.01-0.02 eV. As a consequence, 
at each snapshot there is a substantial amount of pathways over which the CT might occur 
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FIG. 8: Effect of electronic couplings on CT in (Gg)4 (Na+): Average transmission and current 
calculated i) for the full time series of couplings obtained from MD (full), ii) only intra- and inter- 
strand couplings are non-zero (intra+inter) , iii) the quadruplex is separated into its 4 single strands 
(S4s), for which the intrastrand transport is calculated independently and added up afterwards, 
and iv) for Gs multiplied by 4 (4xpG). 

through the quadruplex. These findings are supported by the I-V characteristics in Fig. |8j 
The current for the intra+inter model almost matches that of the full MD with 2.0 and 
2.3 nA, respectively. If we switch off the interstrand couplings in the quadruplex (S4s) the 
current will drop down to 0.08 nA which is more than order of magnitude smaller. A slightly 
larger current of 0.11 nA is obtained for 4xpG. 

For further insights into the different CT properties between G4 and dsDNA we make use of 
conformation analysis. For that purpose, we investigate the amount of conformations which 
dominate the average CT. Additionally, the distribution of transmissions for the multitude 
of conformations is explored for both (Gs)4 and Gs- The results are presented in Fig. [9j 
Panel a) evidently indicates that there are substantially more conformations contributing to 
the average CT in G4 than in poly(G). Consequently, virtually every 10th G4 conformation 
is CT-active (about 3000 out of 30,000), whereas only 128 (again out of 30,000) single 
non-equilibrium structures characterize the average CT in poly(G). 

Second, as demonstrated in the transmission probability distribution functions (PDFs) in 
Fig. , the distribution width for the G4 quadruplex is considerably narrower compared to 
4xpG and E4s. Note the x-axis is scaled logarithmic. In addition, the G4 PDF is significantly 
shifted to higher transmission. This underscores that the majority of poly(G) structures is 
not transmissive, yet a few single conformations are responsible for the average CT. On 
the other hand, for E4s the major part of conformations reveal higher transmission than in 
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FIG. 9: Conformation analysis for (Gg)4 (Na + ) and Gs, a): number of conformations that make up 
90% of the average transmission maximum {Tmax), b): probability distribution functions (PDFs) 
of Tmax for (Gs)4 (Na + ) and the two models 4xpG and S4s as used in Fig. m 



4xpG, but the single dominating conformations are missing. Clearly, this explains that these 
single high-transmissive poly(G) conformations are the reason for a better conductivity in 
4xpG compared to S4s. The higher average Tl values in dsDNA are due to few highly 
conducting conformations. The smaller average of the G4 Tl couplings therefore resembles 
a more stable structure, however, not leading to a higher conductance, as could be argued 
beforehand. The advantage of G4 over dsDNA is due to the existence of non-negligible 
interstrand couplings in G4. The amount of high-transmissive structures is remarkably 
increased compared to double-stranded poly(G) DNA. 

a. Contact Effects In most of the conductance experiments for DNA the molecules are 
only connected with one strand to the respective left and right contacts. For instance, in 
a very recent experiment by Scheer et al. a single stranded nucleotide with sequence 5'- 
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(T*Gs[TTAGGG]3T*)-3' (T* denotes modified thymine residues), which is known to form 
stable quadruplexes, has been attached between two contacts [B7]. Therefore, the question 
arises whether the CT in an all-parallel stranded quadruplex differs if only 2 or even 1 strand 
of the quadruplex are coupled to the left and right contacts, respectively. For that purpose, 
the CT in (Gs)4 (Na+) is calculated for various contact models: i) all 4 strands are connected 
to the left and right electrodes, respectively (1-4), ii) only one strand is connected to both 
contacts (1-intra) and iii-v) one strand is attached to the left electrode while one of the 



remaining is contacted to the right one (12-inter, etc.). As can be seen in Fig. 10, reducing 
the number of connected strands from 4 to 1 leads to an decrease in (T(E)) for the relevant 
energy range of about 1 order of magnitude. However, the transmission for the 1-intra and 
inter models is very similar indicating a minor significance for which strand or strands are 
connected to the contacts. This finding is also supported by the I-V characteristics, for the 
current ranges around 0.2 nA for these models which is roughly one order of magnitude 
smaller as though all 4 strands are attached (1-4 model). 

Interestingly, there seems to be an increase in the transmission fluctuations if the quadruplex 
is contacted through only one strand at each end. Notwithstanding, the transmission plateau 
for the 1-stranded (intra and inter) contacted models is still about 1.5 orders of magnitude 
larger compared to those for poly(G) resulting in an average current which is again 1 order of 
magnitude smaller with 0.025 nA. Thus our results suggest that independent on the various 
contact linking schemes of G4 and dsDNA one might expect a higher conductivity for the 
quadruplex. Certainly, the optimal conductance for all-parallel stranded G4-DNA is ensured 
if all four strands are coupled to the contacts. 

b. Effect of DNA Environment on Transport As discussed above, the major part of the 
dynamical disorder is induced by the QM/MM environment, i.e. the last term in Eq. [7]which 
is built of the MM charges of DNA backbone, solvent and counterions. Previous results have 
indicated that the disorder due to the DNA environment might not only suppress CT in 
homogeneous sequences like poly(G), rather it is able to enhance CT in random sequences 
like the heterogeneous Dickerson dodecamer [20]. Very recent experiments by Scheer et al. 
confirmed this notion, since they found the current for a 31 base pair sequence [68J to be 2 
orders of magnitude smaller in vacuo than in aqueous solution [13]. Therefore, it might be 
interesting to investigate, what happens if we switch of the environment for CT in G4-DNA, 
also with respect to the effect in dsDNA. As can be seen in Fig. [lit a), the transmission 
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FIG. 10: Effect of contacts on CT in (G§)4 (Na + ): Average transmission and current calculated for 
various contact models i) all 4 strands are connected to the left and right electrodes, respectively 
(1-4) ii) only one strand is connected to both contacts (1-intra), iii-v) one strand is calculated to 
the left electrode while one of the remaining is contacted to the right one (12-inter, etc.), and vi) 
comparison with poly(G). 

maximum for (Gs)4 in vacuo is about 2 orders of magnitude larger than with the QM/MM 
environment. Besides, its broadening and also the fluctuations are significantly reduced and 
the plateau is shifted to higher energies due to the neglect of the electrostatic interaction 
with the sodium ions within the quadruplex. The transmission for poly(G) in vacuo shows 
basically the same features, for the broadening is likewise strongly reduced and the maximum 
is located in the same energy range as for G4 in vacuo, although it is 2 order of magnitude 
smaller. In general, the transmission for G4 including the QM/MM environment clearly 
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FIG. 11: Influence of MM environment on CT: Comparison between (Gg)4 (Na + ) and Gs- a) 
average transmission and current, b) number of conformations that make up 90% of the average 
transmission maximum {Tmax) and c) probability distribution functions (PDFs) of Tmax- 

reveals the largest broadening which might indicate that the central sodium ions have an 
additional strong influence on the dynamic disorder due to longitudinal mobility within the 
quadruplex, which is not the case in dsDNA. Despite all that, the current at high voltages 
is larger for G4 with QM/MM than for poly(G) in vacuo. 

Basically, for heterogeneous sequences (not shown here) a reduced transmission in vacuo is 
found, which is caused by energy gaps between A and G states. Notwithstanding, for DNA 
molecules with no static energy gaps like double-stranded poly(G) and G4-DNA (both with 
uniform DNA bases), the QM/MM environment is most likely to increase the dynamical 
disorder, thus will suppress CT compared to the vacuo model. As a result, there is no 
significant difference in the effect of the DNA environment on CT for G4 and poly(G) 



DNA. This is also underscored by conformation analysis given in Panel b) and c) in Fig. 11 



which indicate that there are considerably more CT-active conformations in vacuo than with 
QM/MM environment for both DNA species G4 and poly(G). Moreover, for G4 in vacuo 
nearly every conformation appears to be high-transmissive, since the average maximum 
transmission (Tmax) increases almost linearly with the number of conformations. Thus, the 
CT in vacuo is only marginally affected by single non-equilibrium conformations, rather the 
whole ensemble of G4 conformations seems to be CT active which is demonstrated in the 
PDF of transmission. 
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4- Effect of Channel ions on Transport 



One last issue remains when considering CT in G4-DNA that is the effect of ions within the 
quadruplex channel. The structural influence of these ions has already been addressed in 



detail in Sec. Ill A As could be seen in Sec. SII.B in Ref. [72] we did not find the central ions 
to contribute states in the relevant energy range for hole transfer in G4, therefore, the effect 



will only be investigated electrostatically. In Fig. [12] the average transmission and current 
is shown for the quadruplex simulations of (Gg)4 in absence and presence of either lithium, 
sodium and potassium ions. Obviously, the transmission maximum is only marginally af- 
fected by the presence of different types of central ions. Furthermore, as expected the 
transmission function is shifted to lower energies if central ions are present. However, for 
potassium the transmission is to be found slightly reduced, also the broadening is not as 
large as for the other species. This is also reflected in the PDF of transmission maxima 
which can be found in Fig. S4 in Ref. [72]. 

As a consequence, the average current for the simulation with central potassium ions is half 
as large as for those with lithium and sodium ions which might be attributed to the different 
mobilities of Li + and Na + compared to K + . Interestingly, there is no significant difference 



for CT in absence and presence of ions, although it is known from Sec. HI A that the G4 
molecules without central ions exhibit significant destabilization. Once more this supports 
the notion that the enhanced conductance in G4 may not exclusively be explained in terms 
of higher structural rigidity, i.e. less dynamical disorder, rather it is the multitude of CT 
pathways via interstrand couplings that brings on an increased number of high-transmissive 
conformations. Apparently, those interstrand and in-plane couplings T2 and T3 are not 
altered by the presence of central ions, i.e. by a more rigid quadruplex structure. 



IV. DISCUSSION AND CONCLUSION 



In this work, we have investigated the conductivity of G4 with respect to dsDNA using classi- 
cal MD simulations, combined QM/MM methods to compute CT parameters and Landauer 
theory to compute transmission and I-U characteristics. These properties have been evalu- 
ated for phase-space ensembles of G4 and dsDNA structures. 

The approach adopted here to calculate the transport characteristics at a given time (Lan- 
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FIG. 12: Influence of central ions on CT: Calculation of the average transmission and current for 
(Gg)4 in absence and presence of monovalent central ions Li + , Na + and K + . 

dauer theory) has clear limitations related to the underlying assumption of tunneling trans- 
port. Though this may be an efficient pathway over very short segments of the molecules 
under study, its validity becomes questionable with increasing length. This becomes quite 
evident when investigating sequences as used in recent experiments e.g. by Scheer and 
co-workers |13j . where our methodology would predict currents order of magnitude smaller 
than those found in the experiments. We are currently developing a methodology in order 
to describe transport in complex materials in a more general way [T91 152] . 
Nevertheless, the Landauer calculations reported here give interesting insight into particular 
properties of G4, especially when compared to dsDNA. First of all, G4 with central ions is 
structurally much more stable than dsDNA and has more 7r-contacts along the chain [32J. At 
first sight, one may argue that this leads to improved electrical conduction, however, it turns 
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out that the large fluctuations found for dsDNA lead to highly conducting structures, which 
dominate the transport. Therefore, a dsDNA molecule conducts better than one strand of 
G4, or equivalently, the higher conduction of G4 is not due an increased structural stability 
of its single strands. The phase space of the single strands still contains a vast majority of 
conformations, which are not CT active. It is the ability of G4 to allow for a large number of 
conformations due to the interstrand couplings T2 and T3 which makes this species better 
conducting. If the pathway along one strand is blocked, e.g. when one T^- along the chain 
vanishes, many other conduction channels may be viable due to interstrand hopping. At the 
end, G4 has a much larger number of CT active conformations than four dsDNA (poly(G)) 
offers. This is the basis of the higher conductivity of G4. This advantage is even maintained, 
when contacted only at one G-site, instead of four sites. 

As a second point, dsDNA or G4 may be exerted to strain due to the contacting procedure. 
Here, clearly the higher stability of G4 may help to maintain a conducting conformation, 
while the conduction in dsDNA may be much more easily disrupted [B"7| . 
In a recent work, we have predicted dsDNA with non-uniform sequences (i.e. not poly(G) 
or poly (A)) to conduct better in solution than in gas phase. This is due to the larger 
structural fluctuations in solvent, which introduce conducting structures [20]. This finding 
has been confirmed by a recent experiment [13]. For homogeneous sequences like poly(G) 
or G4 however, a higher conductance in gas phase should be found [71] . Finally, our results 
suggest that the presence of central metal ions within the quadruplex has only a marginal 
impact on hole transfer in G4, although they are vital for the stability and rigidity of G4- 
DNA. 
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